Production of heavy isotopes in transfer reactions by collisions of 238u_j_238u 
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The dynamics of transfer reactions in colhsions of two very heavy nuclei 238u_|_238y studied 
within the dinuclear system (DNS) model. Collisions of two actinide nuclei form a super heavy 
composite system during a very short time, in which a large number of charge and mass transfers 
may take place. Such reactions have been investigated experimentally as an alternative way for 
the production of heavy and superheavy nuclei. The role of coUision orientation in the production 
cross sections of heavy nucUdes is analyzed systematically. Calculations show that the cross sections 
decrease drastically with increasing the charged numbers of heavy fragments. The transfer mecha- 
nism is favorable to synthesize heavy neutron-rich isotopes, such as nuclei around the subclosure at 
N=162 from No (Z=102) to Db (Z=105). 

PACS number(s): 25.70.Jj, 24.10.-i, 25.60.Pj 



The synthesis of superheavy nuclei (SHN) is motivated 
with respect to searching the "island of stability", which 
is predicted theoretically, and has resulted in a number of 
experimental studies. There are mainly two sorts of reac- 
tion mechanism to produce heavy and superheavy nuclei, 
namely, multi-nucleon transfer reactions in collisions of 
two actinide nuclei [l|, Q and fusion-evaporation reactions 
with a neutron-rich nuclide bombarding a heavy nucleus 
near shell closure, such as the cold fusion reactions at 
GSI (Darmstadt, Germany) [sl and the '^^Ca bombard- 
ing the actinide nuclei at FLNR in Dubna (Russia) 
However, the decay chains of the nuclei formed by the 
hot fusion reactions are all neutron-rich nuclides and do 
not populate presently the known nuclei. Meanwhile, the 
superheavy isotopes synthesized by the cold fusion and 
the "^^Ca induced reactions are all far away from the dou- 
bly magic shell closure beyond ^"^Pb at the position of 
protons Z=114-126 and neutrons N=184. Multi-nucleon 
transfer reactions in collisions of two actinides might be 
used to fill the region and examine the influence of the 
shell effect in the production of heavy isotopes. The 
production of neutron-rich heavy or superheavy nuclei 
in low-energy collisions of actinide nuclei was proposed 
initially by Zagrebaev et al. based on the assumption 
that the shell effects continue to play a significant role in 
multi-nucleon transfer reactions [5|. 

The cross sections of the heavy fragments in strongly 
damped collisions between very heavy nuclei were found 
to decrease very rapidly with increasing the atomic num- 



ber y] . Calculations by Zagrebaev and Greiner with 
a model based on multi-dimensional Langevin equations 
showed that the production of the survived heavy 
fragments with the charged number Z>106 is rare be- 
cause of the very small cross sections at the level of 1 
pb and even below 1 pb. However, neutron-rich isotopes 
of Fm and Md were produced at the larger cross sec- 
tion of 0.1 fih. The evolution of the composite system 
in the damped collisions is mainly infiuenced by the in- 
cident energy and the collision orientation. Recently, the 
time dependent Hartree-Fock (TDHF) approach [7| and 
improved quantum molecular dynamics (ImQMD) model 
were also used to investigate the dynamics in collisions 

of 238U+238U. 

In this work, we use a dinuclear system (DNS) model 
to investigate the dynamics of the damped collisions of 
two very heavy nuclei, in which the nucleon transfer is 
coupled to the relative motion by solving a set of mi- 
croscopically derived master equations by distinguishing 
protons and neutrons 0, [l^ . In order to treat the diffu- 
sion process along proton and neutron degrees of freedom 
in the damped collisions, the distribution probability is 
obtained by solving a set of master equations numerically 
in the potential energy surface of the DNS. The time evo- 
lution of the distribution probability P{Zi, Ni, Ei,t) for 
fragment 1 with proton number Zi and neutron number 
A^i and with excitation energy Ei is described by the 
following master equations. 



[dz„N,P{Zi,N[,E[,t) - dz,,N[P{Zi,Ni,Ei,t)] ~ [Agf{e{t))+Af,s{emP{Zi,Ni,Ei,t){l) 
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Here the Wzi,Ni;Z[,Ni (W^Zi,Afi;Zi,Af;) is the mean tran- sition probability from the channel {Zi,Ni,Ei) to 
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{Z[ ,Ni,E[) (or {Zi, N,, El) to {Zi,Ni,E[)), and dz, .n, 
denotes the microscopic dimension corresponding to the 
macroscopic state {Zi, Ni, Ei). The sum is taken over 
all possible proton and neutron numbers that fragment 
Z[,N[ may take, but only one nucleon transfer is con- 
sidered in the model with the relation Z[ = Zi ± 1 and 
N[ = iVi ± 1. The excitation energy Ei is determined by 
the dissipation energy from the relative motion and the 
potential energy surface of the DNS. The motion of nucle- 
ons in the interacting potential is governed by the single- 
particle Hamiltonian 0, • The evolution of the DNS 
along the variable R leads to the quasifission of the DNS. 
The quasifission rate Ag^ and the fission rate A of the 
heavy fragment are estimated with the one-dimensional 
Kramers formula. The local temperature is given by the 
Fermi-gas expression Q = ^ e* ja corresponding to the 
local excitation energy e* and the level density parame- 
ter a = A/12 MeV^ 0. 

In the relaxation process of the relative motion, the 
DNS will be excited by the dissipation of the relative 
kinetic energy. The local excitation energy is determined 
by the excitation energy of the composite system and 
the potential energy surface of the DNS. The potential 
energy surface (PES) of the DNS is given by 



f/({a}) =B(Zi,Ari) 
+y({a}) 



B{Z^,N^)-{B{Z,N)^V^,^{J)\ 

(2) 



with Zi + Z2 = Z and iVi + iV2 = N .11]. Here 
the symbol {a} denotes the sign of the quantities 
Zi,Ari,Z2,7V2; J,i?;/?i,/32,0i,^2. The B{Z,,NO{i = 
1.2) and B{Z,N) are the negative binding energies of 
the fragment {Zi, Ni) and the compound nucleus (Z, N), 
respectively, which are calculated from the liquid drop 
model, in which the shell and the pairing corrections are 
included reasonably. The is the rotation energy of 

the compound nucleus. The f3i represent the quadrupole 
deformations of the two fragments. The 6i denote the 
angles between the collision orientations and the symme- 
try axes of deformed nuclei. The interaction potential 
between fragment {Zi,Ni) and (Z2, iV2) includes the nu- 
clear, Coulomb and centrifugal parts; the detailed cal- 
culations are given in Ref. [10]. In the calculation, the 
distance R between the centers of the two fragments is 
chosen to be the value at the touching configuration, in 
which the DNS is assumed to be formed. So the PES 
depends on the proton and neutron numbers of the frag- 
ments. Shown in Fig. 1 is the calculated PES as func- 
tions of mass numbers of fragments for the two cases of 
the nose-nose and side-side orientations. Dissipation to 
heavier fragment by nucleon transfer is hindered in the 
nose-nose collisions and a pocket appears around the nu- 
cleus '^^^120. The side-side orientation is easily to reach 
the subclosure ^^^No, but is also hindered if further dis- 
sipating the high-mass region. 

The cross sections of the primary fragments (Zi, Ni) 
after the DNS reaches the relaxation balance are calcu- 



lated as follows: 



apr{Zi,Ni) 
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Jmax 

E 

J=0 



(2J + l)P(Zi,iVi,r,„t)- (3) 



The interaction time Tint in the dissipation process of 
two colliding partners is dependent on the incident energy 
Ec.m. in the center-of-mass (cm.) frame, and the angular 
momentum J, which is calculated by using the deflection 



function method [l^ and has the value of few 10 ^"s 



The survived fragments are the decay products of the 
primary fragments after emitting the particles and 7 rays 
in competition with fission. The cross sections of the 
survived fragments are given by 



^ sur 



2 J,r 



^(2J+l)P(Zi,iVi,£;i,T„t) 



J=0 

xWsur{Ei,xn,J), 



(4) 



where the Ei is the excitation energy of the fragment (Zi, 
Ni). The maximal angular momentum is taken as Jmax — 
200 that includes all partial waves in which the transfer 
reactions may take place. The survival probability W^ur 
of each fragment can be estimated by using the statistical 
approach |10| . 

Dynamics of the damped collisions was investigated 
by Zagrebaev and Greiner in detail with a model based 
on multi-dimensional Langevin equations 0, [13] • Larger 
cross sections in the production of neutron-rich heavy 
isotopes in collisions of two actinides were pointed out. 
Within the framework of the DNS model, we calculated 
the production of the survived fragments in collisions of 
238u_|_238u ^jjg nose-nose and side-side orientations 
at 900 MeV center-of-mass energy as shown in Fig. 2. 
The collisions of the side-side case need to overcome the 
higher barrier of the interaction potential in the forma- 
tion of the DNS. But it is favorable to transfer nucleon by 
the master equations in the driving potential and form 
the target-like fragments. Situation is opposite for the 
nose-nose collisions. So both cases give a similar result in 
the production of the Bk, Cf, Es and Fm isotopes. Com- 
parison of the calculated mass distributions and the ex- 
perimental data of the survived fragments at Ec.m. =800 
MeV is shown in Fig. 3. The cross sections decrease 
drastically with the atomic numbers of fragments. The 
calculated results are the case of the nose-nose collisions, 
which have the height of the interaction potential at the 
touching configuration with the value 713 MeV. In the 
collisions of such heavy systems, there is no Coulomb 
barrier in the approaching process of two colliding part- 
ners. A number of nucleon transfers take place in the 
reactions of two actinides owing to the dynamical defor- 
mations and the fluctuations of all collective degrees of 
freedoms in the model of Zagrebaev and Greiner [6| . We 
assumed the DNS is formed at the touching configuration 
in the collisions of two very heavy nuclei. The nucleon 
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FIG. 1. Driving potentials of the tip-tip and waist-waist collisions in the reaction U+ U. 
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FIG. 2. Comparison of the production cross sections of the Bk, Cf, Es and Fm isotopes in collisions of ^^^U with ^ssy ^^j. ^j^g 
tip-tip and waist-waist orientations, respectively. 



transfer is governed by the driving potential in compe- 
tition with the quasifission of the DNS because of the 
collision dynamics. The larger quasifission rate of such 
systems results in the DNS quickly decays into two frag- 
ments because there is no potential pockets. The inner 
excitation energy of the DNS is dissipated from the ki- 
netic energy of the relative motion overcoming the height 
of the interaction potential of two colliding nuclei at the 
touching configuration. The value of the side-side orien- 
tation is 813.5 MeV. Inclusion of all orientations in the 
low-energy damped collisions is important for correctly 
estimating the cross sections of the primary and survived 
fragments. 



Shown in Fig. 4 is the calculated cross sections of the 
primary and survived fragments by using the Eqs (3) 
and (4) as functions of the charged numbers and mass 
numbers at the incident energy Ec.m.=800 MeV, respec- 
tively. In the damped collisions, the primary fragments 
result from a number of nucleon transfer in the relaxation 
process of the colliding partners. The giant composite 
system retains a very short time of several tens 10~^^ 
s due to the strong Coulomb repulsion. Calculations 
from the TDHF method showed that the collision time 
depended on the orientation of the colliding system 0]. 
The cross sections in the production of heavy target-like 
fragments (Z>92) decrease drastically with the atomic 
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of heavy isotopes is discussed systematically. The low- 
energy transfer reactions in the damped collisions of the 
actinide nuclei are a good mechanism to produce the 
neutron-rich heavy isotopes, in which the shell closure 
plays an important role in the estimation of the cross 
sections. 
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FIG. 3. Calculated mass distributions of the Cf, Es, Fm and 
Md isotopes at Ec.m.=800 MeV and compared with the avail- 
able experimental data 01 . 



numbers of the fragments. Therefore, the mechanism of 
the low-energy transfer reactions in collisions of two very 
heavy nuclei is not suitable to synthesize superheavy nu- 
clei (Z>106) because of the smaller cross sections at the 
level of 1 pb and even below 1 pb. Similar results were 
also obtained in Ref. However, the production of the 
survived fragments around the subclosure N==162 has a 
larger cross section. Calculated cross sections as func- 
tions of the mass numbers appear a bump near the iso- 
topes of the subclosure. Experimental works for studying 
the influence of the shell closure in the production of the 
neutron-rich isotopes should be performed in the near fu- 
ture. It is also a good technique to fill the gap of the new 
isotopes between the cold fusion and the ^*Ca induced 
reactions. 

In summary, the dynamics in collisions of the very 
heavy system 238u_|_238y jg investigated within the 
framework of the DNS model. The influence of the 
collision orientations on the production cross sections 
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FIG. 4. Cross sections as functions of the charged and mass numbers of the primary and survived fragments at Ec.m,.=800 
MeV, respectively. 



